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We introduce a novel method for computing the (/x, r)-dependent pressure in continuum QCD; 
and therefrom obtain a complex phase diagram and predictions for thermal properties of the dressed- 
quark component of the system, providing the in-medium behaviour of the related trace anomaly, 
speed of sound, latent heat and heat capacity. 
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I. INTRODUCTION 

Strong interactions in the Standard Model are de¬ 
scribed by quantum chromodynamics (QCD), which is 
supposed to describe a vast array of phenomena, from 
gluon and quark interactions at the highest energies 
achievable with the large hadron collider to the nature of 
nuclear material at the core of a compact star. This last 
challenge initiated the quest to uncover the equation of 
state (EoS) for superdense nuclear matter [Q. The inter¬ 
vening years have seen remarkable activity, highlighted 
by the discovery of a new state of matter, viz. a strongly- 
coupled quark-gluon plasma (sQGP) [^. These efforts 
have delivered a sketch of the QCD EoS in the plane 
spanned by baryon chemical potential (rb) and temper¬ 
ature (T) |, |. 

In vacuum, i.e. in the neighbourhood rb — 0 ~ T, 
QCD exists in a phase characterised by two emergent 
phenomena: confinement and dynamical chiral symmetry 
breaking (DCSB). Confinement is most simply defined 
empirically: those degrees-of-freedom used in defining 
the QCD Lagrangian gluons and quarks do not exist 
as asymptotic states. The forces responsible for confine¬ 
ment appear to generate more than 98% of the mass of 
visible matter [^. This phenomenon is known as DCSB. 
It is a quantum field theoretical effect, which is ex¬ 
pressed and explained via the appearance of momentum- 
dependent mass-functions for quarks and gluons 

[|lOHl^ even in the absence of any Higgs-like mechanism. 

On the other hand, in medium, i.e. as rb and/or T 
are increased beyond certain critical values, the prop¬ 
erty of asymptotic freedom suggests that QCD 

undergoes phase transitions. In the new phases, DCSB 
disappears and/or gluons and quarks are deconfined. In¬ 
deed, the possibility that the transitions are related and 
coincident in the (/x^,F)-plane is much discussed. 

The QCD EoS and related thermal properties are im¬ 
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portant for numerous reasons j^, e.g. they are crucial 
inputs to the hydrodynamic simulations used to connect 
basic theoretical predictions with modern experimental 
data. In this context, the current sketches are not ade¬ 
quate: a better picture is needed in order to understand 
the data and hence the qualities of a sQGP. 

The T % 0 properties of QCD have been scrutinised 
via simulations of the lattice-regularised theory (IQCD). 
This is apparent, e.g. in Refs. which also high¬ 

light problems impeding the extension of lattice meth¬ 
ods to /xb % 0. At this stage, a picture of the phase 
diagram in the entire (/x^, r)-plane requires other meth¬ 
ods, so models continue to be employed widely. However, 
they are numerous in number and various in formulation, 
and too often provide conflicting predictions [|^ . Herein, 
therefore, we analyse aspects of the thermal properties of 
QCD using methods of continuum quantum field theory; 
namely, the Dyson-Schwinger equations (DSEs) 

II. DRESSED-QUARK PRESSURE 

We focus on the matter sector of QCD and hence begin 
with the quark gap equation: 

S{p, 0 Jn)~^ = i^ ■ p -\- i^A^n A rn P E{p,ujn), (la) 
°° f d^a 

^{p,(An) = T I —^g^D^^{p-q,nrii;T,R) 

l = — oo'^ ^ ^ 

X Y7/i*S'(g,a;/Yr^(5,(I;i,:P,a;„), (lb) 

where m is the current-quark bare-mass; P ip, 

LUn = (2n+l)7rr is the quark Matsubara frequency and r 
is the quark chemical potential {r^ = 3/x), Clni = — 

is the dressed-gluon propagator; and Fj^ the dressed- 
quark-gluon vertex. 

The kernel of Eq. (j^) is determined by the quark-gluon 
vertex and the gluon propagator. We use Fj, = ju, which 
defines rainbow-ladder (RL) truncation, i.e. the leading- 
order in the most widely use d sy mmetry-preserving DSE 
approximation scheme |^, |^, which is accurate for 
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ground-state light-quark hadrons It is appro¬ 

priate here because whilst a certain class of vertex im¬ 
provements can influence critical exponents associated 
with second order transitions, viz. those leading to in¬ 
clusion of long-range colour-singlet correlations Q, no 
form available today alters the order of a transition or 
has a material impact on its location p9|^^. 

The gluon propagator in Eq. ( 0 ) has the general form 

= pj:,DT(k\u^) + . (2) 

where are, respectively, k transverse and longi¬ 
tudinal projection operators, and -H P^^, = — 

Dt = P(so,0), Dl = P(so,mp, where = 
(16/5)(r^-|-6/i^/[57r^]) describes a gluon screening mass, 
the value of which is determined from leading-order per¬ 
turbative QCD Q; and (sq = + k‘^ + m?g) Q 


P(5o,mJ) 


7mP(go) 

ln[r+(l + V^QCD)^] 


( 3 ) 


with zF(z) = (1 - = 0.5GeV, r = - 1 , 

7 ^ = 12/25, Aqcd = 0.234GeV. P(sQ,m^) is shape- 
consistent with solutions of in-vacuum gap equations [|4| ; 
and with Du = (0.8 GeV)^ and a renormalisation-group- 
invariant current-quark mass rhu,d = 6 MeV, the solu¬ 
tions of Eq. (|^) support a reliable in-vacuum description 
of ground-state hadrons in RL truncation Q. {N.B. By 
using Landau gauge we minimise sensitvity to truncation- 
induced violations of gauge covariance [Q.) 

The in-medium extension of the gap equation’s ker¬ 
nel, Eq. (||), preserves agreement with QGD at large 
momenta. However, in assuming that D is {ld,T)- 
independent it overlooks screening of the interaction’s 
infrared strength. We remedy that by writing [|^ 

d(t,h) = dIi --. r > r / ( 4 ) 

[ b{^) -\-ln[f /Aqcd] 

where + 6 /i^/[ 57 r^] and Pp marks the onset of 

thermal screening. With Pp = Pc, the critical temper¬ 
ature for chiral symmetry restoration, then a = 0.029, 
b = 0.47 yield a dressed-quark thermal mass ttit = 0.8 P 
at P = 2Pc, in agreement with IQGD simulations [^. 
Naturally, Pc = Tc{fd): we set Pp(/i) = Tc{fj) at p, 7 ^ 0, 
ensuring D{Tc{fd), fi) = D by evolving the value of 5(/i). 

In RL truncation and stationary phase approximation, 
the dressed-quark pressure density is |^ : 

P[S] = Pin Z = -P(TrIn [TS] + ^Prp*?]), (5) 

where Eq. determines S, E. At each {fd,T), Eq. ® 
possesses the same ultraviolet divergence, which may be 
eliminated by subtracting the /r = 0 = P result. The 
subtraction can be accomplished by recalling that for any 



FIG. 1. Points - Numerical result for the free-quark pressure 
obtained as described in connection with Eq. (|); and solid 
(red) curve - analytical result. 


function f{w), compatible with a physical system, Q: 


ruo 

27TiT f{iu^Pii)= / dwf(w 

n= — oo '^“0 


n= —00 
rUiJ.+V 


/u*-\-r] 


dw f{w) 

Q{w-fj,)/T I 




dw f{w) 

u*—ri ^ m J- 


. ( 6 ) 


where we have omitted a P-independent term that is zero 
so long as /i < /ic(T), and = iA -\- fi with A ^ 00 , 
7^0. The hrst term on the right-hand-side of Eq. (^ is 
responsible for the divergence we wish to eliminate. The 
physical piece of the dressed-quark pressure can thus be 
calculated as the difference between the two terms in the 
hrst line of Eq. (^, with f{w) computed via the func¬ 
tional expression in Eq. (||). Practically, one proceeds as 
follows; solve the gap equation at a given (/r, P)-pair for a 
large number of Matsubara frequencies, characterised by 
rim', at each /r, obtain smooth interpolations in “m” for 
the scalar functions obtained thereby, Eq. (|^; evaluate 
P[S] with these inputs; and then compute the difference 
between the sum and integral, verifying that it is insen¬ 
sitive to the interpolation procedure and choice of rim- 
That this procedure can work effectively is illustrated in 
Eig. ^ which displays a comparison of our numerical re¬ 
sult for the free-quark pressure with the analytic form. 

Our approximation for the dressed-quark contribution 
to the QGD EoS is now dehned. Using the pressure, one 
can dehne the trace anomaly: X = e — 3P, where the 
energy density £ = —PpTs and s(T) = dP[S]/dT is the 
entropy density. Notably, P = 0 for an noninteracting ul- 
trarelativistic gas, which is described by Psb oc and 
hence P is a measure of the interaction energy stored in 
the system. Now, since the conhned dressed-quark con¬ 
tribution to the physical pressure must vanish on P ~ 0, 
it follows that P exhibits a maximum at some Tm, the 
value of which serves to dehne a useful reference temper¬ 
ature. 
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FIG. 2. Trace anomaly (normalised, = 0): I/Psb- Psb is 
the pressure of a noninteracting ultrarelativistic gas contain¬ 
ing the number of gluons and quarks appropriate to the calcu¬ 
lation. Solid (black) curve - dressed-quark contribution, com¬ 
puted via Eq. (^. For comparison, IQCD results obtained us¬ 
ing various discretisation schemes. Complete pressure: short- 
dashed Iblue) long-dashed (green) and dot-dashed 
(red) p|]. Gluon-only contribution: dotted (pink) curve [pif . 


III. PHASE DIAGRAM AND 
THERMODYNAMICS 

We employ the DSEs because they possess the capac¬ 
ity to study confinement and DCSB simultaneously in 
the continuum Q. Within this framework the (/x 7 ^ 
0 , T 7 ^ 0) EoS has only been computed using a very sim¬ 
ple description of QCD’s gauge sector |^. Owing to the 
importance of the EoS in developing a complete picture 
of the Standard Model, it is imperative to do better. 

We depict the /x = 0 trace anomaly in Eig. ^ In order 
to facilitate comparisons between the profiles obtained in 
different analyses, the temperature is expressed in units 
of the appropriate value of Tm- The value of Tm ex¬ 
hibits modest variation between the calculations: in our 
case Tm = 0.14 GeV; which is roughly 30% smaller than 
found in modern IQCD studies owing to our omission of 
the gluon contribution. Additionally, each computation 
in Eig.^ is normalised by the appropriate form of Psb: 
represented in this way, there is qualitative agreement 
between all results. Our prediction describes the quark- 
only contribution to X. It matches the IQCD results in 
shape and order of magnitude. These observations indi¬ 
cate both that the gluon and quark contributions to the 
total interaction energy behave similarly and that they 
are commensurate in size when measured against their 
respective asymptotic contributions to the pressure. 

The “speed of sound” in the system is obtained from 
Cg = dPjde. It is a crucial factor in determining the flow 
of material in the system, i.e. its transport properties. 
As evident in Eig.^, our prediction for the sound velocity 
in the dressed-quark subcomponent is similar to results 
obtained in IQCD for the velocity in the complete system. 

A feature which distinguishes our framework from 
IQCD is its ability to treat /x > 0 without further approxi¬ 
mation. Thus, in Eig. § we display results for the dressed- 
quark pressure and trace anomaly at a range of values of 


FIG. 3. Temperature dependence of c^, where Cs is the sound 
velocity in the system. Solid (black) curve - our result, com¬ 
puted from the dressed-quark pressure; points (blue) - IQGD 
results from Ref. 0; and band (green) - IQGD results from 

Ref. Q. 

/X > 0: whilst increasing /x produces an increase in the 
pressure at all values of temperature, it only materially 
increases the interaction energy on T < Tm- Similar be¬ 
haviour is seen in IQCD estimates for the /x-dependence 
of these quantities, obtained using various extrapolation 
schemes ||^, . In detail, P is a monotonic function 

of T for small fi, and P and X are smooth; but qualitative 
changes occur at /x = /Xp = 0.106 GeV: a peak appears in 
P and that in X becomes sharper. Both functions remain 
smooth, however, until /x = 0.111 GeV, whereat the first 
derivative of each diverges at P = 0.128 GeV, signalling 
that the transition has become first-order. This effect lo¬ 
cates the critical endpoint (CEP) for the chiral symmetry 
restoring transition at (/x^ = 0 . 111 , = 0.128) GeV. 

The behaviour of both P and P on fi/T 2 =; 0 is consistent 
with hard thermal loop perturbation theory j^. 

We now draw the diagrams associated with QCD’s 
phase transitions as determined from the dressed-quark 
pressure. Chiral symmetry restoration is straightfor¬ 
ward. It may be charted via the (/x, r)-dependence of 
the chiral condensate but we prefer a method 

Q based on the chiral susceptibility, x(M)T). 

Before discussing deconfinement, however, one must 
have a definition of colour confinement. We consider 
confinement as a violation of reflection positivity by 
coloured Schwinger functions. This associates it with 
dynamically-driven changes in the analytic structure of 
QCD’s propagators and vertices that occur be¬ 

cause both gluons and quarks acquire running mass dis¬ 
tributions, which are large at infrared momenta. This 
leads to the emergence of a length-scale ~ 0.5 fm, whose 
existence and size is evident in all continuum and IQCD 
studies of dressed-gluons and -quarks: characterises the 
material change in their analytic structure 
Erom this perspective, deconfinement occurs when ^ 0 
and reflection positivity is thus recovered. This crite¬ 
rion has been used effectively in dense-hot QCD (e.g. 
Refs. 1^, 1^, and we employ it herein, following 

a local implementation elucidated in Refs. [^, |^, |^ . 

We display the phase diagram computed from the 







4 





FIG. 4. Upper panels: Temperature dependence of the 
dressed-quark pressure {left) and trace anomaly {right) at sev¬ 
eral chemical potentials. Lower left panel: latent heat of the 
chiral symmetry restoring transition. Lower right panel: heat 
capacity of the system’s dressed-quark subcomponent. 



FIG. 5. Deconfinement and chiral symmetry restoration 
phase boundaries, computed via pressure in Eq. (|^. Gurves: 
solid (black) - chiral transition, with DGSB favoured below 
the curve; dashed (green) - deconfinement transition, with 
dressed-quarks confined below the curve; dot-dashed (blue) 
- Nambu phase chiral susceptibility, xn'- positive below the 
curve and zero above; dotted (red) - Wigner phase chiral sus¬ 
ceptibility, xw- positive above the curve and negative below. 
GEPs: chiral - filled circle; and deconfinement - asterisk. 


dressed-quark pressure in Fig.|[ Comparison with Fig. 3 
in Ref. shows that our improved DSE kernel, which 
agrees with the one-loop QCD renormalisation group, 
does not qualitatively alter the phase diagram. The solid 
curve in Fig. ^ is the locus of transition; the Nambu phase 
is energetically favoured for those values of (/x, T) that lie 
within the domain bounded by the axes and this curve; 


Tf(z = = r 


p _ _p 1 + 0.52^2 - 0.058z'‘ 


cO' 


1 + 0.91^2 


( 7 ) 


The /X = 0 pseudo-critical temperature associated with 
the chiral transition is = 0.15 GeV. For comparison, 
a IQCD estimate of the critical temperature for chiral 
symmetry restoration in QCD with two light flavours and 
a physical strange quark mass is Tc = 0.15 ± 0.01 GeV 
[^. Within a factor of two, the z = 0 slope of Tf {z)/T^ 
from Eq. (^ agrees with estimates from IQCD |4l| . 

With (/X, T) increasing from the origin, the dot-dashed 
curve in Fig.|| bounds the domain of positive Nambu- 
phase chiral susceptibility. The dotted curve, on the 
other hand, marks the line whereat the Wigner-phase 
chiral susceptibility switches from negative to positive. 
These curves coincide with the transition locus, Eq. (|^), 
up to the chiral transition’s CEP; 


CEP^ = {fil = 0 . 111 ,TJ = 0.128) GeV, ( 8 ) 

which confirms the result obtained in connection with 
Fig. i top panels; and the computed rati o /x^ 

0.87 is commensurate with those in Refs. 

The chiral crossover becomes a first-order transition at 
CEP;^; the Nambu and Wigner phases coexist, with the 
Nambu phase dominant below the transition locus and 
the Wigner phase dominant otherwise. 

The dashed curve marks the boundary for the decon¬ 



finement transition, which is second-order until 

CEP^ = (/x|; = 0.185,r| = 0.106) GeV, (9) 

/x^/r|, = 1.75. Wigner phase-domains exhibit neither 
confinement nor DCSB at any values of (/x, T). In the 
chiral limit the solid and dashed curves coincide and 
CEP^ = CEP^; chiral symmetry restoration and decon¬ 
finement are coincident; but the dislocation at m 7 ^ 0 
entails the existence of a small domain wherein quarks 
are deconfined but chiral symmetry is broken. This is 
the set {/X, T} enclosed within the dot-dashed and dashed 
curves. Pockets of Nambu phase material in this subset 
of the phase coexistence region possess that character. In 
the domain of (/x,r) enclosed between the axes and the 
dashed and dotted curves, the system is in a pure phase 
with confinement and DCSB, whereas chiral symmetry is 
restored and quarks are deconfined in the domain above 
the dot-dashed curve. 

We can now compute the heat-capacity density of the 
dressed-quark system, cy = dejdT^ and the latent- 
heat density of transition; L = T/S.s = Ae — /xAp, 
p{T) = dP/dp is the quark number density and Ap is 
the difference between the quantity F in the two distinct 
phases, which are here the Nambu (chirally asymmet¬ 
ric) and Wigner (chiral symmetry restored) phases. We 
depict L, computed along the phase boundary, which is 
the trajectory in Eq. (^), in the lower-left panel of Fig.||. 
Naturally, L = 0 for T > Tj because the transition is 
no longer first order. Otherwise our prediction is qual¬ 
itatively consistent in shape with what may be inferred 
from IQCD simulations Idl] . 

The lower-right panel of Fig.^ displays our predic¬ 
tion for cy; it diverges as (/x, T) ^CEP^^. Actually, 
in the neighbourhood of the CEP one may write [^; 
cv oc \g — where g = fi, T; and the quark number 
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susceptibility x = dp/dfj, behaves in the same fashion. 
Analysing our results, we find e = 0.67 ± 0.02 in both 
cases. These are the critical exponents of a mean-field 
transition, which is the nature of RL truncation. 

IV. SUMMARY 

We introduced a practical procedure for computing 
the (/i, T)-dependent dressed-quark pressure in contin¬ 
uum QCD, which we illustrated using a gap equation 
whose solutions are key to a successful description of the 
properties of ground-state hadrons in-vacuum. Without 
further approximation, the associated richly-structured 
phase diagram in the (/r, r)-plane was computed. We 
drew the transition lines for deconfinement and chiral 
symmetry restoration, and confirmed that these transi¬ 
tion are identical in the chiral limit. Likewise, we calcu¬ 
lated the speed of sound in the system and provided the 
(/i, r)-dependence of the trace anomaly, latent-heat and 
heat-capacity densities. Where comparisons are possible, 
our predictions are consistent with results from lattice 
QCD. Notably, predictions obtained from the dressed- 


quark pressure are qualitatively equivalent to those com¬ 
puted using the complete pressure, which suggests that 
the dressed-quark pressure alone can be used as a practi¬ 
cal guide to some of QCD’s thermal properties. No ma¬ 
terial improvement over our results can be envisaged in a 
continuum analysis before a symmetry-preserving kernel 
including long range correlations is derived for the gap 
equation. Our method for computing the pressure will 
also be applicable then. 
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